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ABSTRACT 

A deployment is the movement of armed forces from their 
home bases to their strategic locations. The movement of 
these forces usually involves the transportation of 
military personnel as well as equipment and supplies. In a 
crisis situation, it is essential that the deployment is 
carried out in a expeditious manner. 

This study considers the problem of constructing a 
deployment plan for sealift assets which transport military 
personnel, equipment, and supplies to their designated 
locations in the least amount of time. In the construction 
of such a plan, feasible transportation schedules for each 
asset must be specified. When the number of movement 
requirements is large, the problem of arranging schedules 
for the assets is nontrivial. This thesis, therefore, 
describes an algorithm to generate these schedules. Based 
on several examples, this algorithm is shown to be 
effective and can be used in conjunction with algorithms 
for solving the overall deployment problem as well. 
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I. INTRODUCTION 



A. PROBLEM STATEMENT 

A deployment plan is a collection of schedules, one for 
each of the lift assets (trucks, trains, planes, and ships) 
assigned to the deployment. Each schedule normally 
comprises a list of cargoes arranged in the sequence which 
they are to be picked up from the ports of embarkation 
(POE's) and delivered to ports of debarkation (POD's) by 
the designated asset. 

The generation of schedules becomes a great concern in 
the construction of optimal deployment plans because of the 
enormous number of possible schedules. Even for a 
deployment which consists of only 10 lift assets and 20 
shiploads of cargoes, the number of possible schedules 
exceeds 10 million. Thus, regardless of the type of 
algorithm employed for the overall problem, the ability to 
produce an optimal deployment plan in a timely manner 
greatly depends on the ability to produce a good set of 
feasible schedules efficiently. 

B. BACKGROUND 

In mathematical programming, the idea of generating 
columns for the constraint matrix was first introduced by 
Dantzig and Wolfe [Ref. 1]. Among the first to employ the 
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column (schedule) generation (or the Dantzig-Wolfe 
decomposition) technique in the area of deployment are Rao 
and Zionts [Ref. 2] and Appelgren [Ref. 3]. Rao and Zionts 
employ the technique to solve the problem of allocating 
ships to existing schedules, and Applegren applies it to 
solve a ship scheduling problem. 

More recently, Brown, Graves, and Ronen [Ref. 4] and 
Brown, Goodman, and Wood [Ref. 5] use a slightly different 
approach in solving the problem of scheduling oil tankers 
and naval surface combatants. Instead of generating the 
schedules as needed in the manner of the Dantzig-Wolfe 
decomposition, they have to generate all feasible schedules 
apriori because the schedules must satisfy many rules, 
restrictions, and relationships among the ports, ships, and 
cargoes. However, these rules, restrictions, and 
relationships also reduce the number of acceptable 
schedules to a manageable amount, thereby making the 
generation of schedules feasible. 

C. OBJECTIVE 

Most authors (see, Ronen [Ref. 6], Ronen [Ref. 7], 
Brown, Graves, and Ronen [Ref. 4], Brown, Goodman, and Wood 
[Ref. 5], Collier [Ref. 8], and Lally [Ref. 9]) consider a 
deployment with the minimum cost to be optimal. Instead of 
cost, one can also minimize other cost related quantities 
such as total distance traveled or the number of lift 
assets required. These quantities, if necessary, can always 
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be converted into a monetary amount. However, this study 
adopts a different optimality criterion which arises in the 
military deployment planning during the period of conflict 
(crisis) . In this situation, the speed at which the armed 
forces, equipment, and supplies are deployed to the area of 
conflict is one of the most important factors governing how 
the conflict is resolved in the end. Therefore, it is 
essential that the plan can deliver all cargoes to their 
destinations in the shortest amount of time. 

The objective of this study is to develop an efficient 
algorithm to generate schedules that can be used in 
conjunction with any column generation scheme for solving 
the crisis deployment problem. One such scheme is presented 
in a related Naval Postgraduate School Masters Thesis by 
Lt. N. R. Lima [Ref. 10]. 

The remainder of this thesis is divided into 4 
chapters. Chapter 2 gives a summary of the mathematical 
formulation of the deployment planning problem. Chapter 3 
describes the schedule generator. Chapter 4 discusses the 
computer implementation and the results from the 
experimentation with a group medium size problems. Finally, 
Chapter 5 summarizes the study and indicates areas for 
future research. 
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II. PROBLEM FORMULATION 



In the deployment problem addressed below, all cargoes 
are assumed to be in full shiploads. This implies that the 
asset assigned to pick up a given cargo must deliver it 
before any other cargo can be picked up. Because of this 
assumption, the cargoes are also referred to as "movement 
requirements". Moreover, it is also assumed that other 
data essential to the problem, such as the distances 
between ports, ships' speeds, ships' initial positions, and 
tables indicating the ship/port and ship/cargo 
compatibility, are also given. Then, the problem of 
constructing a deployment plan with the minimum duration of 
time can be stated as follows: 

Indices : 

i - movement requirements (cargoes), where i = 1,....,M 
and M is the number of movement requirements. 

j - ships, where j = 1, ,N and N is the number of 

available lift assets. 

k - feasible ship schedules. As mentioned in Section C of 
Chapter 1, the generation of these schedules is the 
topic of this thesis. 
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Data : 



= 1 if and only if the kth feasible schedule for the 

jth ship includes the ith movement requirement. 
Thus, Sjjr is an M-dimensional vector representing 
the kth feasible schedule for ship j. A schedule 
for ship j is considered feasible if ship j is 
compatible with the cargoes it must pick up and 
with the ports it must visit. 

tjj^ = the completion time of schedule Sj)^. 

Decision variable : 

Xj)^ = 1 if and only if the kth schedule for ship j is 

selected for the deployment. 

Problem PI : 

Min { max ( Etjj^Xjj^ : j = 1,...,N ) } 

j k 

Subject to : 

N 



s 

j=l 


^ ^ijk ^jk - 

k 


for i = 1 , . . . ,M 


(1) 


2 

k 


Xjk ^ 


for j = 1 , . . . ,N 


(2) 




Xjk = 0 or 1 


• 


(3) 


In the 


above formulation 


constraint (1) forces 


each 



movement requirement to be picked up by at least one ship, 
while constraint (2) allows each ship to sail at most one 
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schedule. A set of x's satisfying constraints (1), (2), and 

(3), if it exists, would specify a deployment plan to 
support all the requirements of the deployment using the 
given assets. 

In a related thesis, Lima [Ref. 10] obtains a solution 
to Problem PI by parametrically solving a simpler problem 
several times. This simpler problem. Problem P2 (t) below, 
is called the feasibility-seeking problem and is a variant 
of Problem PI. In the feasibility-seeking problem, one 
assumes that a maximum time for deployment, t, is given. 
Then, the task is to find a set of schedules with 
completion time less than or equal to r that satisfies the 
constraints on Problem PI. To state the problem 
mathematically, let 

Auxiliary Variables : 

Wj^ = 1 if and only if movement requirement i is not 

delivered, and 

Auxiliary Data : 

Kj(T) = the index set of feasible schedules for ship j 
with the completion time less than or equal to 

T , i. e. , 

Sj(T) = { k : Sj}^ is a schedule and tjj^ < t }. 

Then, the feasibility-seeking problem can be stated as 
follows: 
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Problem P2 (t) : 



Min E Wj^ 
i 

subject to : 
N 



^ ^ ^ijk ^jk + > 

j=l keKj(T) 


1, for i 


= 1,...,M 


(4) 


S Xjk ^ 

k 


1, for j 


= 1,...,N 


(5) 


Xj]^ = 0 or 


1. 




(6) 


Associated with each constraint 


in (4) 


there is a 


dual 



variable (<0) which is referred to as the "cargo duals" 
in the subsequent chapters. Similarly, associated with each 
constraint in (5) there is a dual variable Vj (>0) which 
will be referred to as the "ship dual". 

Observe that for a given value t, if the optimal 
objective function value of P2 (t) is zero, then the upper 
bound for Problem PI is t. On the other hand, if the 
optimal objective function value of P2 (t) is positive, then 
there is at least one undelivered shipload of cargo. 
Therefore, t must be increased in order to obtain a 
feasible plan with the same number of assets. By varying t 
and resolving Problem P2(r) in a systematic manner, one can 
obtain a solution to PI. 

Although Kj(r) is a subset of the feasible schedules, 
it is still too large to enumerate. However, one can avoid 
this enumeration of schedules in the set Kj(r) by applying 
the Dantzig-Wolfe decomposition technique to Problem P2 (t) 
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which produces a master and a subproblem. The master 
problem has the same form as Problem P2 (t) and the 
subproblem generates schedules (or columns) which belong to 
the set Kj ( T ) and have negative reduced costs . 

The following chapter describes an efficient algorithm 
for the subproblem. The details of solution techniques and 
strategies for the master problem are presented in the 
related thesis by Lima [Ref. 10] . 



8 



III. A SCHEDULE GENERATOR 



A. THE SUBPROBLEM 

The subproblein resulting from the application of the 
Dantzig-Wolfe decomposition to Problem P2 (t) can be stated 
as follows: 

Problem Sl(j) : 

k(j) = arg min { vj + 2 : tjj^ < r } 
k i 

As previously defined in Chapter 2, u^ and Vj are the cargo 
and ship duals of Problem P2 (t). For a given deployment 
plan, a negative cargo dual, Uj^, indicates that the ith 
shipload is undelivered, and its magnitude can be regarded 
as the profit obtained from delivering the cargo to its 
destination. The ship dual, Vj , can be interpreted as the 
cost of using ship j. Therefore, the problem Sl(j) would 
find a new schedule which utilizes ship j in the most 
profitable manner. 

Moreover, Problem SI is indexed by j to indicate that 
the generated schedule is for ship j . In theory one can 
solve one subproblem for a particular ship j , or solve N 
subproblems, one for each ship. To insure convergence only 
one schedule (column) with a negative reduced cost needs to 
be added to the master problem during each cycle of the 
decomposition process. 



9 



If the set of schedules Sjj^ is available apriori, 
solving Problem Sl(j) is merely a matter of selecting one 
schedule from the set Kj(r). Otherwise, Problem Sl(j) 
resembles the vehicle routing problem [Ref. 11] with the 
exception that, 

(1) the ship schedules do not have to start and 
terminate at the same port, 

(2) the shiploads of cargoes must be delivered before 
any other shiploads can be picked up, and 

(3) the schedule must be completed before time t. 

These three additional conditions greatly increase the 
complexity of Problem Sl(j). 

As stated. Problem Sl(j) would always produce a 
schedule with the most negative reduced cost for each ship 
j. Moreover, considering the complexity of the problem, it 
may not be advantageous to solve Problem Sl(j) to 
optimality. Moreover, the convergence of Dantzig-Wolfe 
decomposition only requires that the reduced cost of the 
new schedule (column) be negative. Thus, one only needs to 



search for a schedule. 


Sjjr, which satisfies the 


following 


two conditions: 






Vj + S 

i 


^ijk ^i <0 


(7) 




tjk ^ ^ 


(8) 



where condition (7) ensures that the schedule has a 
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negative reduced cost and (8) guarantees that the schedule 
can be completed within the allowable time, t. 

In practice, entering columns which merely satisfy 
conditions (7) and (8) does not generally lead to an 
efficient algorithm. In the early iterations of the 
Dantzig-Wolfe algorithm, the ship duals, v j , tend to be 
zero and the cargo duals are mostly negative. Thus, a 
schedule which picks up only one cargo with a negative dual 
would easily satisfy both conditions (7) and (8) . However, 
such a schedule would be considered a bad schedule if the 
ratio of cargoes (movement requirements) to the number of 
ships is much larger than 1. For example, if the ratio is 
four then on the average one would expect each ship to 
complete four movement requirements. So, a schedule that 
picks up only one cargo is unlikely to be included in the 
optimal plan. The algorithm described below will search for 
one or more schedules which pick up a minimum number of 
cargoes as well as satisfying conditions (7) and (8) . This 
minimum number of cargoes is calculated based on the ratio 
of cargoes to ships. 

B. AN ALGORITHM WITH THE DEPTH-FIRST SEARCH STRATEGY 

A schedule which satisfies conditions (7) and (8) and, 
at the same time, picks up the minimum number of cargoes is 
considered to be a "good" schedule. The algorithm below 
searches for one or more good schedules in the same manner 
as the branch and bound algorithm, with the depth-first 



11 



strategy, searches for an optimal solution for an integer 
program. 

Basically, the algorithm tries to generate a sequence 
of numbers which is called a path. With the exception of 
the number zero, which represents the ship's initial 
position, the numbers in the path represent the cargo 
number and the order of the numbers represents the sequence 
in which the cargoes are to be picked up and delivered. 
Thus, a schedule or column for the master problem can 
easily be constructed from this path. 

To illustrate the principle underlying the algorithm, 
consider a deployment problem with 3 movement requirements 
numbered 1 to 3 . Figure 3-1 depicts a tree which represents 
all possible schedules for a ship. The node number in the 
tree corresponds to the movement requirement number except 
for the root node, which is numbered zero and corresponds 
to the ship's initial position. By traversing this tree 
starting at the root node, one would trace out a path from 
which a schedule can be constructed. For example, consider 
path 1 in Figure 3-1. This path corresponds to the 
following sequence of numbers: 0-1-3 -2 which in turn 
corresponds to a schedule for the ship to complete 
movement requirements 1,3, and 2 in sequence. Similarly, 
path 2 corresponds to a schedule for the ship to complete 
movement requirements 2 and 3 in sequence. 
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Initially, the algorithm discards all movement 
requirements which are either incompatible with the ship in 
consideration or have a positive dual variable. This 
incompatibility is due to either ship/port or ship/cargo 
incompatibility. Assume that the remaining movement 
requirements are numbered as 4 , 6, 7, and 9. Then, starting 
at node zero with the current path containing only the 
number 0, the algorithm creates nodes 4, 6, 7, and 9 and 

attaches them to node zero (see Figure 3-2) . From node 
zero, the algorithm considers adding a new movement 
requirement to the current path, say node 4; this is called 
branching. At node 4, the algorithm checks to see if the 
schedule which contains movement requirement 4 can be 
completed in t days. If the answer is yes, the algorithm 
then includes node 4 into the current path, creates nodes 
corresponding to movement requirements which are not 
members of the current path, and attaches them to node 4. 
These nodes are node 6, 7, and 9. At this point, the 

algorithm would branch to, say, node 6 and at node 6 it 
checks to see if the schedule which contains movement 
requirements 4 and 6 can be completed in r days. If no, the 
node is not added to the current path. The algorithm would 
then check if the current path, 0-4, corresponds to a 
schedule with negative reduced cost and picks up the 
minimum number of shiploads of cargo. If no, node 6 is 
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"fathomed" and the algorithm would branch to other 
unfathomed nodes in the depth-first manner, which in this 
case would be node 7 or node 9 in the second level. If yes, 
then a good schedule is found, and the algorithm would stop 
if only one good schedule is desired. If more schedules are 
required, the algorithm would store this good schedule, 
fathom the node, and branch to other unfathomed nodes to 
obtain the desired number of schedules or until all nodes 
in the tree are fathomed. 

Below, we formally state the algorithm discussed above. 
It employs the abstract data type called STACKS [Ref. 12] 
to dynamically generate the search tree and conserve memory 
requirements . 



Algorithm 



In put ; 

j = the index of the ship being considered. 

T = the maximum allowable completion time for any 

schedule. 



M = 

NCOL = 
NPICK = 




the total number of shiploads of cargo 
(movement requirements) . 

the required number of "good" schedules for 
ship j . 

the minimum number of cargoes to be picked up 
by each good schedule. 

the dual price for ship j . 

the dual price for movement requirement i . 
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Cj^j = 1 if movement requirement i is compatible 

with ship j 

= 0 otherwise. 

Step 1 : For each i = set 

Ui = Ui*Cij + 99*(1 - Ci j ) . 

Remove all movement requirements i with 
Ui>0 from further consideration, and let M* be 
the remaining movement requirements, 
i. e., those with u^sO. 



Step 2 : Sort in an ascending order. 

Step 3 : Set KOUNT = 0, P = {0}, and initialize the 

STACK. 



Step 4 : 



Step 5 : 



5.1) 



5.2) 



Put the movement requirements into the 
STACK in the same order as the sorted dual 
variables u^. 

While the STACK is not empty, do the 
following: 

Remove the movement requirement, say, k, 
the top of the STACK. 

Calculate the completion time for path 



from 



P U {k}. 

5.3) If this completion time s g, then do the 
following: 
a) Set P = P U {k}. 
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b ) Put the movement requirements not in the path 



on top of the STACK in the order of decreasing 

Ui- 

c ) Go to Step 5 . 1 

5.4) Otherwise, calculate the reduced cost for the 
schedule corresponding to path P. 

a) If the reduced cost is negative and the 
number of movement requirements in P r NPICK, 
then a good schedule is found. Set KOUNT = 
KOUNT + 1. 

If KOUNT s NCOL, remove node k from the STACK 
and go to Step 5.1. 

Otherwise, return to the master problem. 

b) If either the reduced cost is nonnegative or 
the number of movement requirements in P is 
less than NPICK, then remove node k from the 
STACK and go to Step 5.1. 

Note that Step 5.3.b dynamically generates the new 
level of nodes in the search tree. Moreover, the order in 
which the nodes are put on top of the STACK allows the 
search to progress in the depth-first manner with the 
movement requirement with the most negative cargo dual 
being searched first. In Step 5.4, the removal of node k 
from the STACK is equivalent to fathoming the node in the 
search tree. 
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IV. COMPUTATIONAL RESULTS 



The algorithm for generating schedules presented in 
Chapter 3 was implemented in VS FORTRAN on the IBM 3033AP 
computer at the W. R. Church Computing Center of the Naval 
Postgraduate School. The master problem is also implemented 
in the same manner and is fully described in the Master's 
thesis by Lima [Ref. 10]. 

For the computational experiments, we consider three 
problems which contain approximately the same number of 
lift assets and movement requirements as one would expect 
in a deployment of a small armed force such as that of 
Norway. Problem 1 has 48 movement requirements and 25 lift 
assets, Problem 2 has 48 movement requirements and 30 lift 
assets, and Problem 3 has 60 movement requirements and 30 
lift assets. The ports of embarkation (New York, Norfolk, 
Charleston, Jacksonville, Pensacola) and debarkation 
(Hamburg, Wilhelmshaven, Rotterdam, Antwerpen, Cherebourg) 
for the movement requirements are given in Table 4-1. Table 
4-2 and 4-3 give the distances between the various ports 
and the initial position of each asset. The speed for the 
assets are between 13 and 22 knots, and the probability 
that the ship is compatible with a given movement 
requirement is .75. 
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TABLE 4-1 : NUMBER OF MOVEMENT REQUIREMENTS 
BETWEEN POE'S AND POD'S. 





HAM. 


WILH. 


ROT. 


ANTW. 


CHB. 


p 1 
R 

0 a 
B n 
L d 
E 

M 2 


N. Y. 


2 


2 


2 


2 


2 


NORF. 


3 


2 


2 


2 


2 


CHAR. 


2 


2 


2 


2 


2 


JAX. 


2 


3 


3 


2 


3 


PENS. 


0 


0 


2 


2 


0 





HAM. 


WILH. 


ROT. 


ANTW. 


CHB. 


P 

R 

0 

B 3 
L 
E 
M 


N. Y. 


4 


4 


3 


2 


2 


NORF. 


2 


2 


2 


3 


2 


CHAR. 


3 


3 


2 


2 


2 


JAX. 


3 


2 


2 


2 


3 


PENS. 


0 


1 


3 


4 


2 



TABLE 4-2 : DISTANCE BETWEEN POE'S AND POD'S 
IN NAUTICAL MILES. 





HAM. 


WILH. 


ROT. 


ANTW. 


CHB. 


N. Y. 


4030 


3950 


3790 


3775 


3520 


NORF. 


4340 


4260 


4090 


4075 


3800 


CHAR. 


4650 


4560 


4390 


4370 


4090 


JAX. 


4850 


4770 


4590 


4570 


4280 


PENS. 


5390 


5300 


5125 


5110 


4820 
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TABLE 4-3. 



INITIAL DISTANCE BETWEEN SHIPS AND PORTS 



SHIP # 


N. Y. 


NORF. 


CHAR. 


JAX. 


PENS. 


1 


0 


245 


550 


720 


1190 


2 


0 


245 


550 


720 


1190 


3 


300 


545 


850 


1020 


1490 


4 


300 


300 


500 


700 


1100 


5 


800 


800 


900 


1000 


1400 


6 


100 


350 


650 


820 


1290 


7 


245 


0 


300 


475 


975 


6 


245 


0 


300 


475 


975 


9 


200 


200 


400 


600 


1100 


10 


600 


600 


700 


900 


1400 


11 


150 


100 


. 400 


575 


1075 


12 


350 


100 


200 


375 


875 


13 


550 


300 


0 


165 


700 


14 


550 


300 


0 


165 


700 


15 


800 


750 


700 


750 


1100 


16 


450 


200 


100 


265 


900 


17 


350 


250 


300 


350 


750 


18 


1200 


1100 


1000 


1100 


1500 


19 


720 


475 


165 


0 


600 


20 


720 


475 


165 


0 


600 


21 


llOO 


1000 


900 


900 


1200 


22 


920 


675 


365 


200 


400 


23 


450 


350 


300 


400 


800 


24 


1350 


1250 


1200 


1200 


1600 


25 


1190 


975 


700 


600 


0 


26 


1190 


975 


700 


600 


0 


27 


890 


675 


400 


300 


300 


28 


1290 


1075 


800 


700 


300 


29 


900 


700 


500 


400 


600 


30 


1090 


875 


600 


500 


100 
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A. A HEURISTIC FOR SELECTING MOVEMENT REQUIREMENTS 

In Step 1, the algorithm discards those movement 
requirements which have a positive dual value because of 
incompatibility with the ship under consideration. The 
remaining movement requirements are then selected for 
branching based solely on the magnitude of the associated 
cargo duals. Intuitively, this seems inefficient. Assume 
that the algorithm is constructing a new schedule for ship 
j . Then those movement requirements which are in the 
schedule for ship j and those which are still unassigned at 
the end of the preceding master problem should be selected 
for branching first. Otherwise, movement requirements which 
have been assigned to other ships may be assigned to the 
new schedule for ship j also. This would create a "double 
coverage" for the same movement requirement. To create the 
least amount of double coverage a negative number is added 
to the cargo duals for the movement requirements which are 
already assigned to ship j and to those which have yet to 
be assigned to any ship. This negative number should be 
large enough to insure that when the movement requirements 
are ordered in Step 2 and put into the STACK in Step 3 and 
Step 5.3.b, the desired movement requirements are always 
above the others in the STACK. 

To implement the above heuristic. Step 1 of the 
algorithm in Chapter 3 is modified as follows: 
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step 1 : 



a) For each i = 1 ,..., M set 



99 if c^j = 0. 



Uj^ if c^j = 0 and cargo i has been 
assigned to another ship 1 , 

_ where 1 = j , by the previous 

Ui = < master problem iteration. 



6 if Cj^j = 0 and cargo i is 

unassigned or has been assigned 
to ship j by the previous 
master problem iteration. 



where 5 is a sufficiently large positive number, 
b) Remove all movement requirements i with 

Uj^>0 from further consideration, and let M' 
be the remaining movement requirements, i.e., 
those with ui<0. 

The algorithm in Chapter 3 has been implemented both 
with and without the above heuristic for selecting movement 
requirements. Figure 4-4 shows that, based on the average 
CPU time to solve the three sample problems with various 
values for r , the algorithm with the heuristic clearly 
outperforms the one without. In fact, the algorithm with 
the heuristic is on the average 47% faster. Henceforth, the 
subproblem algorithm is implemented with the heuristic for 
selecting movement requirements. 
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PERCENTAGE OE OPTIMAL DEPLOYMENT DURATION 



Figure 4-4 

Comparison of total CPU time with and without 
heuristic implemented. 
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B. THE NUMBER OF SCHEDULES TO GENERATE FOR THE MASTER 

PROBLEM 

In general, it is unclear what is the optimal number of 
columns to generate during each cycle of the decomposition 
process. More columns mean more information for the master 
problem during each cycle. At the same time, more columns 
in the master problem also mean more CPU time for pricing 
out and other housekeeping operations. However, Figure 4-5 
demonstrates that it is advantageous to generate between 10 
and 20 schedules (columns) for Problem 2 when r = 29 days 
(the optimal duration) . 

Table 4-6 provides an explanation for the above 
finding. The table lists the average number of nodes in 
the search tree for the two strategies; one column and ten 
columns. Although the average number of nodes in the search 
tree for the ten column strategy is more than the one 
column strategy in the early iterations, the information 
provided by the additional nine columns actually helps the 
master problem to find a solution quicker. Moreover, as the 
number of iterations grow, the one column strategy 
occasionally generates unacceptably large search trees in 
order to find one negative reduced cost column which could 
have been generated earlier if the ten columns strategy was 
employed. 
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NUMBER OF SCHEDULES GENERATED 



Figure 4-5 

CPU time for subproblem versus 
number of schedules generated. 
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TABLE 4-6 



: AVERAGE NUMBER OF NODES GENERATED BY THE 

ALGORITHM FOR PROBLEM 2. 



One column per cycle 


Ten columns per cycle 


Iter 


Average # 
nodes 


Iter 


Average # 
nodes 


1-10 


143 


1-10 


459 


11-20 


104 


11-20 


402 


21-30 


104 


21-30 


414 


31-40 


178 


31-40 


1482 


41-50 


149 


41-45 


315 


51-60 


139 






61-70 


1323 






71-80 


331 






81-90 


297 






91-100 


4999 






101-110 


152 






111-120 


13889 






121-130 


290 






131-140 


475 






141-150 


19233 






151-160 


8751 






161-169 


483 
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Figure 4-7 and 4-8 depict the trade-off between time 
spent in the master and subproblem for Problem 2. Figure 4- 
7 shows that by generating more columns per iteration one 
uses less CPU time for the subproblem. However, more 
columns per iteration also means more columns for the 
master problem which in turn leads to more CPU time for 
pricing out and other-housekeeping operations. In terms of 
total CPU time for the overall problem, Figure 4-8 shows 
that the 10 column strategy is better than the one column 
strategy when t is less than or equal to the optimal 
duration (29 days) . However, when t is bigger than 29 days, 
the one column strategy is superior. This is due partly to 
the fact that when t is larger than the optimal duration, 
the problem is more relaxed, i.e., it has a larger feasible 
region, and accurate dual information is not as critical as 
in the case when t is smaller than the optimal duration. 

Table 4-9 summarizes the computational results for all 
three problems. It is clear from this table that the ten 
column strategy dominates the one column strategy for all 
three problems. However, it would be premature to make any 
conclusive recommendations based on the results presented 
here. 
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DEPLOYMENT DURATION IN DAYS 



Figure 4-7 

CPU time for subproblem versus t 
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DEPLOYMENT DURATION IN DAYS 



Figure 4-8 

Total CPU time versus t 
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TABLE 4-9 



COMPUTATIONAL RESULTS FOR THE THREE PROBLEMS 





PROBLEM 1 


PROBLEM 2 


PROBLEM 3 


1 col 


10 col 


1 col 


10 col 


1 col 


10 col 


7 


CPU 

time 


CPU 

time 


CPU 

time 


CPU 

time 


CPU 

time 


CPU 

time 


.8*OPT 


307 


94 


69 


42 


413 


406 


•9*OPT 


356 


206 


234 


143 


550 


477 


OPT 


470 


177 


112 


9 


2048 


447 


l.l*OPT 


18 


5 


5 


4 


4 


5 


1.2*OPT 


4 


2 


3 


4 


10 


5 


Average 


231 


97 


85 


40 


605 


268 
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V. CONCLUSIONS AND FUTURE RESEARCH 



An algorithm for generating schedules for the crisis 
deployment problem is described. This algorithm is similar 
to the branch and bound algorithm for integer programs or 
the depth-first search technique for traversing directed 
graphs. The computational results in Chapter 4 demonstrate 
that on medium size problems the algorithm, when integrated 
with the master problem algorithm by Lima 10] » is 

effective in solving the subproblem of the crisis 
deployment problem. However, a more extensive computational 
study is required before any conclusive recommendations can 
be made. 

In addition to further computational studies, other 
possible areas for future research are listed below. 

(1) Allow movement requirements to arrive at POE's within 
a time window. 

(2) Allow movement requirements in partial shiploads. 
This would imply that a ship can pick up cargoes at 
two or more locations before delivering them to their 
discharging ports. 

(3) Investigate other branching strategies for the 
algorithm. For example, one can branch to the closest 
movement requirement instead of the movement 
requirement with the most negative dual value. 

(4) Investigate other combinatorial algorithms for the 
subproblem. 

(5) Consider possible delays due to random events such as 
attacks on ships or ports. 
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APPENDIX 



FORTRAN CODE FOR SUBPROBLEtt 



SUBROUTINE SUBPR(U, XCOL, HDAY, h, N, NLOA, NPOE. NPOD, NSH, HR, TR, IT, A, 
&COMPAT, IB, XB, SHIP, XTIHE, K, SPD, SEO, LSEO, CTSHIP, CTNODE, NPICK, NCOL ) 





THIS SUBROUTINE APPENDS NEW (GOOD) COLUMNS TO THE A MATRIX 


IN 


THE 






MASTERPROBLEM. 










KEY VARIABLES COMING FROM THE MASTER PROBLEMi 










- Ui 


DUAL VARIABLES. 










- MDAYl 


MAXIMUM DURATION OF ANY SCHEDULE. 










- Ai 


A MATRIX. 










- COMPATl 


SHIP COMPATIBILITY WITH SHIPLOADS OF CARGO. 










- IBi 


INDEX SET FOR THE BASIS. 










- XBi 


COLUMNS IN THE BASIS. 










- NPICKi 


MINIMUM NUMBER OF SHIPLOADS OF CARGO TO BE PICKED 


BY 








ANY SHIP. 










- NCOLi 


NUMBER OF SCHEDULES TO GENERATE FOR EVERY CALL 


TO 


THE 








SUBPROBLEM. 










KEY VARIBLES USED IN THE SUBPROBLEMi 










- VINDi 


INDEX SET OF THE SORTED DUAL VARIABLES. 










- LOADi 


SHIPLOAD NUMBER. 










- PATHj 


SEQUENCE OF SHIPLOAD NUMBERS FOR THE CURRENT SCHEDULE 






- STACK 1 


SHIPLOAD NUMBERS IN THE STACK. 










- CURLDi 


CURRENT SHIPLOAD NUMBER. 










- PREDi 


PREDECESSOR OF SHIPLOAD NUMBER. 










- TIMEj 


COMPLETION TIME FOR GENERATED SCHEDULES. 










- LENGTH 1 


LENGT OF CURRENT PATH. 









IMPLICIT DOUBLE PRECISION (A-H,0-Z), INTEGER (I-N) 

PARAMETER ( MM -100, NN - 2000, KK - 2, JJ - 2000 ) 

DIMENSION XCOL(MM), U ( MH ) , UU ( HM ) , V ( MM ) , XB ( MM ) , SA ( MH ) , SPD (MM), IB (MM) 
REAL A(MM, NN) 

INTEGER VIND(MM), PRED(Oi JJ), LOAD(Oi JJ), TIME(Oi JJ), STACK(Oi JJ), TOP, 
i MR( 100, KK), TR( 15, 15), IT(30, 5), FROLD, TOLD, PATH(Oi MM), CURLD, COUNT, 

S. LENGTH, LASTND, SHIP, TT, MLNGTH, XTIME(NN), CTSHIP, CTBACK ( 0 i MM ) , 

& CTSON(OiMM), SEO(NN, MM), LSEQ(NN), CTNODE 
LOGICAL COMPAT(MM, MM) 

DOUBLE PRECISION MIN, MINRC 
COMMON /UNITS/ NIN, NOUT 

C INITIALIZE 

CTNODE - 0 
LIMIT - N 
CTSHIP - 0 

311 DO 360 I = 1, MM 
SA( I ) = 0. ODO 
UU( I ) - 0. DO 
V( I ) « 0. DO 

VIND(I) - 0 

360 CONTINUE 

NNEG • NLOA 
MINRC • O.ODO 



SORT THE DUAL VARIABLES 



DO 8 I - 1, NLOA 
8 UU( I ) - U( I ) 



C HEURISTICi 

C ASSIGN LARGE NEGATIVE DUAL VALUE TO SHIPLOADS NOT ALREADY PICKED 
C BY OTHER SHIPS 



202 

201 



IF( N . GE. (2-NLOA ♦ NSH ♦ 1) ) THEN 

DO 201 I - 1, M 

DO 202 J - 1, M 

IF( (IB(J) .GT. (2-NLOA ♦ NSH)) .AND. 

A (XB(J) .GT. 0.5D0) ) THEN 

SA(I) - SA(I) ♦ A(I, IB(J) ) 

END IF 
CONTINUE 

IF( SA(I) . LT. l.ODO ) UU(I) - -2. ODO 

CONTINUE 

END IF 
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C REMOVE SHIPLOADS NOT COMPATIBLE WITH SHIP FROM FURTHER CONSIDERATION. 
DO 9 I ■ 1, NLOA 

9 IF( .NOT. COMPAT( I. SHIP) ) UU(I) ■ 99. ODO 

DO 20 I ■ 1. NLOA 

MIN ■ 0. lD-6 
COUNT - 0 
IND ■ 0 

DO 30 J - 1, NLOA 

IF( UU(J) . LT. MIN ) THEN 
MIN ■ UU(J) 

IND ■ J 
COUNT ■ 1 
END IF 
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CONTINUE 






IF( COUNT .EQ. 
NNEG ■ I - 
GO TO 40 
END IF 
V(I) ■ MIN 
VIND(I) ■ IND 


0 ) THEN 

1 


20 


UU(IND) » 99. ODO 
CONTINUE 




40 


CONTINUE 





• DEPTH FIRST SEARCH 

• — — — — — — — — — — - — — ^ - — — 

C INITIALIZE 

DO 350 I -0. JJ - 1 
TIME(I) ■ 0 
STACK(I) ■ 0 
LOAD(I) - 0 
PRED(I) ■ 0 
350 CONTINUE 

DO 380 I ■ 1, MM 
XCOL( I ) ■ 0. DO 
PATH(I-l) ■ 0 
CTBACK(I-l) - 0 
CTSON(I-l) ■ 0 
380 CONTINUE 



LENGTH ■ 0 
LASTND ■ 1 
TOP ■ 0 
CURLD - 0 
FROLD • 0 
TOLD - 0 

C CREATE ALL NODES OUT OF THE SOURCE AND PUT THEM IN STACK 

DO 50 I ■ NNEG. 1. -1 

LOAD(NNEG - I ♦ 2) - VIND(I) 

PRED(NNEG - I ♦ 2) ■ 1 

STACK(NNEO - I ♦ 1) - (NNEG - I ♦ 2) 

TOP - TOP ♦ 1 
LASTND ■ LASTND ♦ 1 

50 CONTINUE 

C MAIN LOOP TO SEARCH FOR FEASIBLE SCHEDULES 

100 CURLD - STACK(TOP) 

C SPECIAL CASE WHEN ONE SHIP CAN PICK ALL SHIPLOADS OF CARGO 

IF( LENGTH . EQ. NNEG ) THEN 
RCOST ■ 0.000 
DO 91 I ■ 1. LENGTH 

91 XCOL(LOAD(PATH( I ) ) ) - 1. ODO 

XCOL(NLOA ♦ SHIP) - 1. ODO 
DO 111 I ■ 1, M 

111 IF( XCOL(I) . EQ. l.ODO ) RCOST - RCOST ♦ U(I) 

IF(RCOST .QT. -l.OD-4) RCOST - 0. ODO 
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IF( (RCOST . LT. O.ODO) ) THEN 
N • N^l 

DO 389 I - 1, « 

369 A( I, N) • XCOL( I ) 

DO 122 J-i, LENGTH 

122 SEO(N. J) • LOAD(PATH( J) ) 

LSEO(N) - LENGTH 
RETURN 
END IF 

END IF 

C IF STACK EMPTY RETURN TO MASTERPROBLEM IF GOOD SCHEDULE FOUND. 

IF( CURLD .EG. 0 ) THEN 

IF( ttINRC . LT. O.ODO ) RETURN 

C IF NO GOOD SCHEDULES FOUND RUN SUBPROBLEtt AGAIN FOR ANOTHER SHIP. 

CTSHIP ■ CTSHIP ♦ 1 
SHIP ■ SHIP ♦ 1 

IF<SHIP .EG. NSH ♦ 1) SHIP - 1 
IF( CTSHIP .NE. NSH ) GO TO 311 
RETURN 
END IF 

TOP ■ TOP - 1 

IF( PRED(CURLD) .EG. 1) THEN 
LASTND ■ LASTND - LENGTH 
LENGTH ■ 0 
END IF 



IF( CTSON(LENGTH) .EG. CTBACK ( LENGTH ) ) THEN 

DO 51 I - 1, LENGTH 

IF( PRED(CURLD) .EG. PATH(I) ) THEN 
LASTND •= LASTND - LENGTH I 
LENGTH - I 
GO TO 539 
END IF 

51 CONTINUE 

END IF 

539 PATH(LENGTH ♦ 1) • CURLD 

FROLD - LOAD(PRED(CURLD) ) 

TOLD ■ LOAD(CURLD) 

TT • 0 

C CALCULATE COMPLETION TIME FOR CURRENT SCHEDULE. 

CALL TIMEFC(NPOE, NPOD, MR, TR. IT, FROLD, TOLD, SHIP, NSH, NLOA, TT, SPD) 
TIME(CURLD) - TI ME ( PRED ( CURLD ) ) ♦ TT 

C IF FEASIBLE INCLUDE SHIPLOAD NUMBER IN PATH. 

IF( TIME(CURLD) . LE. MDAY ) THEN 

CTBACK(LENGTH) ■ CTBACK ( LENGTH ) ♦ 1 
LENGTH • LENGTH ♦ 1 
CTSON( LENGTH) ■ 0 

C PUT INTO THE STACK ALL SHIPLOAD NUMBERS NOT ALREADY IN THE PATH. 

DO 60 I • NNEG, 1, -1 

DO 70 J • 1, LENGTH 

70 IF( VIND(I) .EG. LOAD ( PATH ( J ) ) ) GO TO 60 

LASTND ■ LASTND ♦ 1 
LOAD(LASTND) - VIND(I) 

PRED(LASTND) - CURLD 
TOP ■ TOP ♦ 1 

CTSON< LENGTH) ■ CTSON ( LENGTH ) ♦ 1 



60 CONTINUE 

DO 60 I - LASTND, (LASTND - CTSON ( LENGTH ) ♦ D,-l 
STACK(TOP) • I 
TOP - TOP - 1 
60 CONTINUE 
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TOP ■ TOP ♦ CTSON( LENGTH) 



ELSE 

LASTND - LASTND - 1 
RC05T • O.ODO 

C CHECK IP GOOD SCHEDULE FOUND. 

DO 90 I 1, LENGTH 

90 XCOL(LOAD<PATH(I) ) ) ■ 1. ODO 

XCOL<NLOA ♦ SHIP) ■ l.ODO 
DO 110 I - 1, M 

110 IF< XCOL<I) . EO. l.ODO ) RCOST - RCOST ♦ U<I) 

IF<RCOST . GT. -l.OD-4) RCOST ■ O.ODO 
IF< RCOST .LT. O.ODO .AND. LENGTH . GT. NPICK ) THEN 



IF < CTBACK(LENGTH) . EQ. 0 ) THEN 

N - N ♦ 1 
DO 3B8 I - 1, M 

368 A< I, N) - XC0L( I ) 

XTIttE(N) - TIHE(PREDCCURLD) ) 
DO 322 J-1, LENGTH 

322 SEO<N,J) ■ LOAD<PATH( J) ) 

LSEQ(N) - LENGTH 

IF< RCOST .LT. MINRC ) THEN 
MINRC - RCOST 
K - N 
END IF 

END IF 

DO 400 I • 1, LENGTH 

XCOL(LOAD(PATH<I) ) ) - O.ODO 
XCOL<NLOA ♦ SHIP) ■ O.ODO 
400 CONTINUE 

ELSE 

DO 130 I - 1, LENGTH 

130 XCOL<LOAD(PATH< I ) ) ) ■ O.ODO 

XCOL(NLOA ♦ SHIP) - O.ODO 

END IF 

C IF NO GOOD SCHEDULE FOUND BACKTRACK, 

CTBACK<LENGTH) - CTBACK ( LENGTH ) ♦ 1 
GO TO 100 

END IF 

IF< N .GE. LIMIT^NCOL) RETURN 
CTBACK(LENGTH) • 0 
GO TO 100 

END 
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SUBROUTINE TIHEFC(NPOE, NPOD, MR, TR, IT, FROLD, TOLD, SHIP, NSH, NLOA, TT, 
iSPD) 

• — 

• THIS SUBROUTINE CALCULATES THE TIME REQUIRED TO PICK UP AND DELIVER • 

• A SHIPLOAD OF CARGO IN WHOLE DAYS. • 



IMPLICIT DOUBLE PRECISION ( A-H, 0-Z ), INTEGER ( I-N ) 

PARAMETER( KK • 2, MM -100, NN • 2000 ) 

DIMENSION SPD(MM) 

INTEGER MRdOO, KK), TR(15, 15), IT(30, 5), TT, TOLD, SHIP, FROLD 
TT - 0 

C CALCULATING THE TRAVEL TIME 
IF(FROLD . EQ. 0) THEN 

TT ■ IDNINT( < IT(SHIP, MR(TOLD, 1) ) ♦ TR ( MR ( TOLD, 1 ) , MR ( TOLD, 2 ) ) ) / 
& (24. • SPD(SHIP) ) ) 

ELSE 

TT • IDNINT( (TR(MR(TOLD, 1 ), MR(FROLD, 2) ) ♦ TR ( MR ( TOLD, 1 ) , 

I MR(TOLD, 2) ) )/<24. • SPD(SHIP))) 



END IF 

RETURN 

END 
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